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Abstract. A new weak Galerkin (WG) method is introduced and analyzed for the second order 
elliptic equation formulated as a system of two first order linear equations. This method, called 
WG-MFEM, is designed by using discontinuous piecewise polynomials on finite element partitions 
with arbitrary shape of polygons/polyhedra. The WG-MFEM is capable of providing very accurate 
numerical approximations for both the primary and flux variables. Allowing the use of discontinuous 
approximating functions on arbitrary shape of polygons/polyhedra makes the method highly flexible 
in practical computation. Optimal order error estimates in both discrete H 1 and L 2 norms are 
established for the corresponding weak Galerkin mixed finite element solutions. 
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1. Introduction. Weak Galerkin (WG) methods refer to general finite element 
techniques for partial differential equations in which differential operators are approx- 
imated by their weak forms as distributions. In 20], a weak Galerkin method was 
introduced and analyzed for second order elliptic equations based on weak gradients. 
In this paper, we shall develop a new weak Galerkin method for second order ellip- 
tic equations formulated as a system of two first order linear equations. Our model 
problem seeks a flux function q = q(x) and a scalar function u = u(x) defined in an 
open bounded polygonal or polyhedral domain il C K d (d — 2, 3) satisfying 

(f.f) aq + Vu = 0, V-q = /, in Q 

and the following Dirichlet boundary condition 

(f .2) u — —g on dil, 

where a = (ctij(x.))dxd £ {L°°(Cl)} d is a symmetric, uniformly positive definite matrix 
on the domain fl. A weak formulation for (|1 . (|1 .2[) seeks q £ H(div,Q) and u £ 
L 2 (fl) such that 

(1.3) (aq,v)-(V-v,u) = (ff,vn>an, Vv £ H(div,Ct) 

(1.4) (V-q, to) = (/,«;), Vw £ L 2 (tt). 

Here L 2 (f2) is the standard space of square integrable functions on !J, V • v is the 
divergence of vector- valued functions v on H(div, Q) is the Sobolev space consisting 
of vector-valued functions v such that v £ [L 2 (d)] d and V • v £ L 2 (Cl), (•, •) stands 
for the L 2 -inner product in L 2 (fl), and (•, -)qq is the inner product in L 2 (dVt). 
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Galerkin methods based on the weak formulation (jl.3[) - (|1.4[) and finite dimen- 
sional subspaces of H(div, CI) x L 2 (Cl) with piecewise polynomials are known as mixed 
finite element methods (MFEM). MFEMs for (|l.i p -(|1.2 p treat q and u as unknown 
functions and are capable of providing accurate approximations for both unknowns 
[D [3j 03 H H E M \M 02] ■ All the existing MFEMs in literature possess local mass 
conservation that makes MFEM a competitive numerical technique in many applica- 
tions such as oil reservoir and groundwater flow simulation in porous media. On the 
other hand, MFEMs are formulated in subspaces of H(div, CI) x L 2 (Cl) which requires 
a certain continuity of the finite element functions for the flux variable. More pre- 
cisely, the flux functions must be sufficiently continuous so that the usual divergence 
is well-defined in the classical sense in L 2 (Cl). This continuity assumption in turn 
imposes a strong restriction on the structure of the finite element partition and the 
piecewise polynomials defined on them [IB] 13 IB] ■ 

The weak Galerkin method introduced in [2DJ was based on a use of weak gradients 
in the following variational formulation: find u G H (CI) such that u — —g on dCl and 

(1.5) (a- 1 V U ,V0) = (/,0), V^eH^Cl), 

where 1 (SI) is the Sobolev space consisting of functions for which all partial deriva- 
tives up to order one are square integrable, Hq(CI) is the subspace of iJ 1 (SI) consisting 
of functions with vanishing value on dCl. Specifically, the weak Galerkin finite element 
formulation in [2DJ can be obtained from (jl.5|) by simply replacing the gradient V by 
a discrete gradient Vd defined by a distributional formula. The discrete gradient op- 
erator Vd is locally-defined on each element. It has been demonstrated [2"Dl [T51 UM [T7] 
that the weak Galerkin method enjoys an easy-to-implement formulation that is pa- 
rameter free and inherits the physical property of mass conservation locally on each 
element. Furthermore, the weak Galerkin method has the flexibility of using discon- 
tinuous finite element functions, as was commonly employed in discontinuous Galerkin 
and hybridized discontinuous Galerkin methods [2] H3] . 

The goal of this paper is to extend the weak Galerkin method of [3DJ to the vari- 
ational formulation (|1.3[l - (|1.4[) by following the idea of weak gradients. It is clear that 
divergence is the principle differential operator in (|1.3[l - (jl.4[l . Thus, an essential part 
of the extension is the development of a weakly-defined discrete divergence operator, 
denoted by (Vd - )) for a class of vector-valued weak functions in a finite element set- 
ting. Assuming that there is such a discrete divergence operator (Vd - ) defined on a 
finite element space Vh for the flux variable q, then formally one would have a WG 
method for (|1.1[) - (|1.2[) that seeks q^ G Vh and Uh G Wh satisfying 

(aq/,,,v) - (Vd • v,Uh) = (g,v ■ n)en, Vv g V h 
(V d • q h , w) = (/, w), \fw G W h , 

where Wh is a properly defined finite element space for the scalar variable. The rest of 
the paper will provide necessary details for a rigorous interpretation and justification 
of the above formal WG method, which shall be called WG-MFEM methods. 

The WG-MFEM methods have the following features. First of all, the finite ele- 
ment partition of the domain CI is allowed to consist of arbitrary shape of polygons for 
d = 2 and polyhedra for d = 3. Secondly, the flux approximation space Vh consists of 
two components, where the first one is given by piecewise polynomials on each poly- 
gon/polyhedra and the second is given by piecewise polynomials on the edges/faces 
of the polygon/poly hedra. The second component shall be used to approximate the 
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normal component of the flux variable q on each edge/face. Furthermore, the scalar 
approximation space Wh consists of piecewise polynomials on each polygon/polyhedra 
with one degree higher than that of the flux. For example, the lowest order of such 
elements would consist of piecewise constant for flux and its normal component on 
each edge/face plus piecewise linear function for the scalar variable on each poly- 
gon/polyhedra. There is no continuity required for any of the finite element functions 
in WG-MFEM. 

One close relative of the WG-MFEM is the mimetic finite difference (MFD) 
method, see [HUH] an< ^ ^ e reference cited therein. Both WG-MFEM and MFD share 
the same flexibility of using polygonal/polyhedral elements of the domain. MFD ap- 
proximates the flux by using only piecewise constants on each edge/face, and the scalar 
variable by using another piecewise constant function on each polygonal/polyhedral 
element, while WG-MFEM provides a wide class of numerical schemes with arbitrary 
order of polynomials. Allowing arbitrary shape for mesh elements provides a conve- 
nient flexibility in both numerical approximation and mesh generation, especially in 
regions where the domain geometry is complex. Such a flexibility is also very much 
appreciated in adaptive mesh refinement methods. This is a highly desirable feature in 
practice since a single type of mesh technology is too restrictive in resolving complex 
multi-dimensional and multi-scale problems efficiently |14j . 

The paper is organized as follows. In Section [2 we present a discussion of weak 
divergence operator in some weakly-defined spaces. In Section[3j we provide a detailed 
description and assumptions for the WG-MFEM method. In Section SI we prove that 
the WG-MFEM scheme has one and only one solution for the diffusion problem 
(|1.2j) . In Section [SI we define some local projection operators and then derive some 
approximation properties which are useful in error analysis. In Section [Bl we shall 
establish an optimal order error estimate for the WG-MFEM approximations in a 
norm that is related to the L 2 for the flux and H 1 for the scalar function. In Section[7l 
we derive an optimal order error estimate in L 2 for the scalar approximation by using 
a duality argument as was commonly employed in the standard Galerkin finite element 
methods |I21 [5] . Finally, we provide some technical results in the appendix that are 
critical in dealing with finite clement functions on arbitrary polygons/polyhedra. 

2. Weak Divergence. The key in weak Galerkin methods is the use of weak 
derivatives in the place of strong derivatives in the variational form for the underlying 
partial differential equations. For the mixed problem ([Lip with boundary condition 
(|1.2|) . the corresponding variational form is given by (|1.3I) and (11.41) . where divergence 
is the only differential operator involved in the formulation. Thus, understanding 
weak divergence is critically important in the corresponding WG method. The goal 
of this section is to introduce a weak divergence operator and its approximation by 
using piecewise polynomials. 

Let K be any polygonal or polyhedral domain with interior K° and boundary 
dK . A weak vector-valued function on the region K refers to a vector- valued function 
v = {vo, Vf,} such that vo £ [L 2 (K)] d and V(,-n e H^^ (dK), where n is the outward 
normal direction of K on its boundary. The first component vo can be understood as 
the value of v in the interior of K, and the second component v;, represents v on the 
boundary of K. The requirement of v& • n e H~?(dK) indicates that we are merely 
interested in the normal component of v?, . Note that Vf, may not be necessarily related 
to the trace of Vo on dK should a trace be well defined. Denote by V(K) the space 
of weak vector- valued functions on K; i.e., 

(2.1) V(K) = {v = {v , v b } : v e [L 2 (K)] d , v b • n e H~* (dK)}. 
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Following the definition of weak gradient introduced in |20j . we define a weak diver- 
gence operator as follows. 

Definition 2.1. The dual of L 2 (K) can be identified with itself by using the 
standard L 2 inner product as the action of linear functionals. With a similar inter- 
pretation, for any v G V(K), the weak divergence o/v is defined as a linear functional 
Vd • v in the dual space of H l (K) whose action on each tp £ H (K) is given by 

(2.2) (Vd • v, ip) K := -(v , V<p) K + (v b • n, ip)d K , 

where n is the outward normal direction to dK , (vo, V(^)x = J K ^o ■ V^AfC is the 
action of vq on W<p, and (v(, • n,ip)gK is the action o/v;, • n on ip G H^(dK). 

The Sobolev space [i? 1 (K )] d can be embedded into the space V(K) by an inclusion 
map iv ■ —?• V{K) defined as follows 

«v(q) = {qk,q|9K>- 

With the help of the inclusion map iy, the Sobolev space [H 1 (K)] d can be viewed 
as a subspace of V(K) by identifying each q G [H 1 (K)] d with iy(q)- Analogously, a 
weak vector-valued function v = {v , v?,} G V(if) is said to be in [H 1 (K)] d if it can 
be identified with a function q G through the above inclusion map. It is 

not hard to see that Vd • v = V • v if v is a smooth function in [_ff 1 (_fi')] d . 

Next, we introduce a discrete weak divergence operator by approximating (Vd • ) 
in a polynomial subspace of the dual of H 1 ^). To this end, for any non- negative 
integer r > 0, denote by P r (K) the set of polynomials on K with degree no more than 
r. A discrete weak divergence operator, denoted by Vd,r - > is defined as the unique 
polynomial (Vd,r • v) G P r (K) that satisfies the following equation 

(2.3) (Vi, r -v,$ K = -(v o ,V0)K + (v 6 -n,0) aK , V</> G P r (K). 



3. Weak Galerkin MFEM: Assumptions and Algorithm. Let be a 

partition of the domain O consisting of polygons in 2D and polyhedra in 3D. Denote 
by Eh the set of all edges or flat faces in 7ft, and let £® = £h\dft be the set of all 
interior edges or flat faces. For every element T £ Ti, we denote by \T\ the area or 
volume of T and by h^ its diameter. Similarly, we denote by |e| the length or area of 
e and by h e the diameter of edge or flat face e G £h- We also set as usual the mesh 
size of 7ft by 

h = max hr- 

Ten 

All the elements of Tk are assumed to be closed and simply connected polygons 
or polyhedra. We need some shape regularity for the partition 7ft described as below. 

Al: Assume that there exist two positive constants g v and g e such that for every 
element T G 7ft we have 

(3.1) g v h d < \T\, ee h d - X < M 

for all edges or flat faces of T. 
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A2: Assume that there exists a positive constant k such that for every element 
T G Tfi we have 

(3.2) Khr < h e 

for all edges or fiat faces e of T. 

A3: Assume that the mesh edges or faces are flat. We further assume that for every 
T G Th, and for every edge/face e G dT, there exists a pyramid P(e,T,A e ) 
contained in T such that its base is identical with e, its apex is A e G T, and 
its height is proportional to Ht with a proportionality constant a e bounded 
away from a fixed positive number a* from below. In other words, the height 
of the pyramid is given by a e hr such that a e > a* > 0. The pyramid is also 
assumed to stand up above the base e in the sense that the angle between 
the vector x e — A e , for any x e G e, and the outward normal direction of e is 
strictly acute by falling into an interval [0,0q] with 6q < |. 

A4: Assume that each T G Th has a circumscribed simplex S(T) that is shape 
regular and has a diameter hs(T) proportional to the diameter of T; i.e., 
hs(T) < 7* hr with a constant 7* independent of T. Furthermore, assume 
that each circumscribed simplex S(T) interests with only a fixed and small 
number of such simplices for all other elements T <ETh- 

Reader are encouraged to check out [10] for a similar, but different type of shape 
regularity assumption for the underlying finite element partition of the domain. The 
shape regularity is required for deriving error estimates for locally defined projection 
operators to be detailed in later sections. 

Let us introduce a space of weak vector-valued functions on Th as follows 

(3.3) V := {v = {v , v b }, v| T G V(T), T E T h } , 

where v|t := {v |t, ^bldr} is the restriction of v on T. In other words, V consists of 
functions defined on elements T GTh and edges/faces e G £h so that its restriction on 
each T and dT is a weak vector- valued function in V(T). Roughly speaking, the space 
V is obtained by patching each local space V(T) continuously through the information 
on dT should such a patching is well defined. 

The weak divergence operator as defined in (|2.2j) can be extended to any weak 
vector-valued function v G V by taking weak divergence locally on each element T . 
More precisely, the weak divergence of any v G V is defined element-by-element as 
follows: 

V d • v = V d ■ (v| r ), on T G %■ 

Similarly, the discrete weak divergence as defined in (|2.3|) can be extended to V by 
defining 

(3.4) V d>r • v = V d , r • (v| T ), on T G Th- 

The definition of weak divergence of v = {vo, v;,} G V requires the value of v on 
each element T, namely vo, and the normal component of v& on each edge or face 
e G £h- Thus, it is the normal component of v on each e£4 that really enters into 
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the equation of discussion in numerical methods. For convenience, we introduce a set 
of normal directions on Eh as follows 

(3.5) T>u = {n e : n e is unit and normal to e, e G 4}> 

In the rest of this paper, we will be concerned with a subspace of V in which the 
second component of v = {vo,Vf,} represents the normal component of v on each 
e £ E h ; i.e., (v b )| e = (v| e • n e )n e . 

A discrete weak vector-valued function v = {vo, Vf,} refers to weak vector-valued 
functions in which both vo and v b are vector-valued polynomials. Since the second 
component is represented as v;, = Vbii e where n e is the prescribed normal direction 
to e 6 4i then vi, is required to be a polynomial on each edge or flat face of dT. 
Recall that, for each element T £ Th, Pk(T) denotes the set of polynomials on T with 
degree no more than k > and Pi(e) is the set of polynomials on e £ Eh with degree 
no more than £ > 0. 

We now introduce two finite element spaces which are necessary for formulating 
our numerical schemes. The first one corresponds to the scalar (or pressure) variable 
defined as follows 

W h = {w £ L 2 (f>); w\ T £ Pk+i(T)}, 

where k > is a non-negative integer. The second one corresponds to vector-valued 
functions and their normal components on the set Eh of edges or flat faces, and is 
given by 

(3.6) V h = {v = {v ,v b } : v \T&[Pk(T)] d ,y b \e=v b n e ,v b eP k (e), e £ £ h } . 

The pair Vh X Wh forms a finite element approximation space for the unknowns q 
and u of the problem (|1.3|) - (I1.4I) . For simplicity of notation and discussion, we shall 
refer the above defined finite element spaces as ([Pk(T)] d , P k (e), Pk+i(Tj) element. 
The lowest order of such element makes use of piecewise constant for the flux variable 
on each element T and its edges/faces and piecewise linear for the pressure (scalar) 
variable. 

The discrete weak divergence (Vd,fe+i*) as defined in (|3.4p and (|2.3j) then provides 
a linear map from the finite element space Vh to Wh- In particular, for any v £ Vh 
and w £ Wh, we have the following relation 

(Vd,fe+i • v, w) : = ^ (Y<j,fc+i • v, w) T 
Ten 

(3.7) = ( I ( v b ■ n)wds - f vo- (Vw)dT j . 

With an abuse of notation, we shall use (V^-) to denote the discrete weak divergence 
operator (Vd,fc+i') m the rest of this paper. 

With the discrete divergence given by (|3.7|) . one might naively formulate a finite 
element method by using (|1.3j) and (|1.4p as follows. Find q^ = {qo,qh} £ Vh and 
Uh E Wh such that 



(3.8) (aqo,v ) - (V«j • v,u h ) = (g,v b ■ n) 9n , Vv = {v , v b } £ V h 

(3.9) {V d -q h ,w) = (f,w), VweWh- 
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Unfortunately, due to an insufficient enforcement on the component q^, the resulting 
system of linear equations from (|3.8[) - (|3.9[) generally does not have a unique solution. 
One remedy to this problem is to stabilize the bilinear form (aqo,vo) by requiring 
some communication between qo and q&. To this end, we introduce three bilinear 
forms as follows. 

(3.10) a (77, v) = (a?7o, v ), r/,veV h , 

(3.11) b(ri,w) = (V d -rj, to), r) eV h ,w eW h , 

(3.12) s(rj, v) = p E /it((?7o - »7&) • n, (v - v b ) ■ n) ax , 77, v G Vh, 

Ter h 

and a stabilized bilinear form a s (-, •): 

a s (*7,v) := 0(17, v) + s(r?,v). 

Here p > is any parameter and hr is the size of T. In practical computation, one 
might chose p = 1 and substitute /it by the mesh size ft. for quasi-uniform partitions; 
i.e., partitions for which hx/h is bounded from below and above uniformly in T. 

Weak Galerkin MFEM Algorithm 1. Let % be a shape regular finite ele- 
ment partition of Q and Vh x Wh fee £/ie corresponding finite element spaces consisting 
of ([P k {T)} d ,P k (e),P k+1 (T)) elements with k > 0. An approximation for |Tn ]) -fO j) 
is given by seeking = {qo,qfj} € and G smc/i £/ia£ 

(3.13) a s (q h ,v) - 6(v, tt ft ) = (g, v 6 • n) a n, 
(3-14) 6(q^, «;) = (/, to), 

/or any v = {vo,Vf,} G and to G WV XTie pair 0/ solutions (q^;tt^) is called a 
weak Galerkin mixed finite element approximation of il. 

4. WG-MFEM: Solution Existence and Uniqueness. For any interior edge 
or flat face e £ £h, let 2\ and Ti be two elements sharing e in common. We define 
the jump of to 6 Wh by 



to 



w \ aTl - w \ot 2 i e e£°, 

to, e G <9£1 



Let us introduce two norms that can be used to measure the accuracy of the solution 
q h and u h obtained from the WG-MFEM ([333]) - (pTH)! . The first norm is defined by 

(4.1) |||v||| 2 := a(v, v) + s(v, v), v G V 

and the second one is defined by 

(4-2) HliV= E HVHII+E^H^HIIe, 

for any to G Wh + H 1 (fl) , where Qb is the standard L? projection operator from L 2 {e) 
into Pfc(e). It is not hard to see that || ■ ||i^ defines a norm in Wh since ||to||i,/j = 
would lead to to = const on each element T and [to] = Qtjto] = on each edge or 
flat face e?4- It follows that to = on each element T. As to § • |||, assume that 
HI v I) = for some v G Vh] i.e., 

(av , v ) + p E ^t((vo - Vft) • n, (v - v 6 ) • n) aT = 0. 

TeT h 
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It follows that v = on each element T and (v — v/,) • n = on each edge or flat 
face e G £h- Thus, we obtain v;, • n = on each e G £h- Since is a vector that is 
parallel to n on each e G £h, then we must have = on each e G £h- The other 
properties for a norm can be easily verified for || • ||i^ and ||| • |||. 

Lemma 4.1. The weak Galerkin mixed finite element method (WG-MFEM) 
H3.13\) - (3.14\ ) has a unique solution. 

Proof. It suffices to show that zero is the only solution of (|3. 13ft - (|3. 14ft if / = g = 0. 
To this end, let / = g = and take v = q and w — Uh in ()3.13[) - (|3.14[) . By adding 
the two resulting equations we obtain 

I q/, |[ = a s (ci h ,q h ) = 0, 

which implies that q/ t = 0. Since g = and q/ t = in (|3.13|) . we have from p.7| that 

(4.3) = -(V d • v, u h ) = (( v °' V ^) T _ < Vb ' n ' Uh ^ dT ^ ' V v e Vh - 

T£T h 

By letting v;, = and vo = Vu^ in (|4.3p . we obtain Vu^ = which implies that Uh 
is a constant on each element. Thus, the equation (|4.3[) leads to 

= ^2 (v b • n, Uk)ar, V v g V/ l7 
which can in turn be rewritten as 

(4.4) 0= ^<|v 6 |,K]) e , VveVi, 

provided that [u^] is dchncd in an order that aligns with the prescribed set of direc- 
tions Vh- Here |vt| stands for the magnitude of the vector v&. By selecting a test 
function v with |v&| = 1 and zero elsewhere we arrive at 

{it h} = on each e G £h- 

This, together with the fact that Uh = const on each element, shows that Uh = on 
the whole domain fl. This completes the proof of the lemma. □ 

5. Local Projections: Definition and Properties. Let us introduce two 
projection operators into the hnitc element spaces Vh and Wh by using local L 2 
projections. For each element T G Th, denote by Qq the projection operator from 
L 2 (T) onto Pk(T). Similarly, for each edge or fiat face eG^, let Qfc be the projection 
operator from L 2 (e) onto Pfc(e). For any v = {vo,Ubn e } G V with vi, G L 2 (e) on each 
edge/face e, denote by Q^v = {Qo( v o), Qb(vb)n e } the L 2 projection of v in Vh- 
Furthermore, denote by Lh the L 2 projection from L 2 (fl) onto Wh- Observe that Lh 
is in fact a composition of locally defined L 2 projections into the polynomial space 
Pfc+i(T) for each element T e Th- 
in the usual finite clement error analysis, one often reduces the error for finite 
element solutions into the error between the exact solution and an appropriately 
defined local projection or interpolation of the solution. For the WG-MFEM method 
discussed in previous sections, this refers to the error between the exact solution and 
its local L 2 projection. The difficulty in estimating the projection error arises from 
the fact that the finite element partition Th contains arbitrary polygons or polyhedra 



9 



that are different from the usual simplices as commonly employed in the standard 
finite element methods |12) . 

For simplicity of notation, we shall use < to denote less than or equal to up to 
a constant independent of the mesh size, variables, or other parameters appearing in 
the inequality. 

Lemma 5.1. Let Th be a finite element partition of Q satisfying the shape regu- 
larity assumptions Al - A 4 as given in Section^ Then, we have 

(5-1) ^||q-Qo4<ft 2!s+1) ||q|| 8 2 +1) 0<s<fc 

TGTh 

(5-2) J2 ll V (<l-3oq)|||</> 2 1|q||* + i, 0<s<k 

Ten 

(5.3) (\\u-L h u\\ 2 T + h 2 T \\V(u-L h u)\\ 2 T ) < h 2( - s+ V\\u\\ s+ i, < s < k + 1. 

Ten 

Proof. To derive (|5.1[) . let S(T) be the circumscribed simplex of T on which q 
can be defined by smooth extension if necessary. Let Qoq be the projection of q in 
the element defined on S(T). It follows that 

(5.4) ||q-Qoq||T < ||q-Qoq||r < ||q - Qoq||s(r) < /4 +1 ||q|U+i,s(T) 
Using the above estimate we obtain 

(5-5) J2 \\<l-Q^\\T<h 2{s+1) £ Nls+W). 

TeT h Ter h 

It follows from the assumption A4 that the set of the circumscribed simplices {S(T) : 
T e Th} has a fixed and small number of overlaps. Thus, the following estimate holds 
true 

X! II<iL+i,s(t) < NU+i- 

Ten 

Substituting the above inequality into (|5.5|) yields the desired estimate (|5.1I) . 

To derive (|5.2[) , we use the triangle inequality and the standard error estimate on 
S(T) to obtain 

||V(q-Q q)||T < ||V(q-Q q)|U ( T) + ||V(Q q - Qo<l)\\s(T) 
(5-6) < Ch s T \\q\\ s+lM T) + Ch^WQoq - Q q||s(r), 

where we have also used the standard inverse inequality in the second line. Notice 
that the assumption A3 on Th implies that there exists a ball B C T with a diameter 
proportional to hx- Thus, we have from Lemma I A. 31 (see Appendix) that 

|]Qoq - Qoq||s(T)< IIQoq - Qoq||s < HQoq - Qoqllr 
< l|q- Qoq||r + ||q- Qoqllr 

<^ +1 ||q|| s+1 , s(T) . 
Substituting the above estimate into (|5.6I) yields 

l|V(q-Q q)||T<^||q|| s+ i,s(T). 
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Summing up the above estimate over T £ 7~h leads to the desired estimate (|5.2[) . 
Finally, the estimate (15.31) can be established analogously to (|5.1I) and (15.21) . □ 

We shall derive two equations that play useful roles in the error analysis for the 
WG-MFEM. The first equation is given by the following Lemma. 

Lemma 5.2. For any q £ let Q^q £ V;, be the projection given by local 

L? projections. Then, on each element T, we have 

(5.7) (V<j ■ (Q/,q), w) T = (V • q, w) T - (q ■ n - Q b (q ■ n), w) 9T 

for all w £ Pk+\(T). Moreover, by summing J<5. 7[ ) over all T £ 7ft we obtain 

(5.8) b(Q h q, w) = (V • q, w) - ^ (q ■ n - Q b (q ■ n), w)g T . 

Ten 

Proof. For q £ [iJ 1 (fi)] d , we apply the definition of Qh and the weak divergence 
V<j ■ (QhQ) to obtain 

(Vd • Qh<l, w) T = -(QoQ, Vw) T + (Qf,(q • n e )n e • n, w) dT 

(5.9) = -(Q q, Vw) T + (Qb(q • n), w) aT 

for all w G Pk+i(T). Here we have used the fact that n e = ±n. Since Qo is the L 2 
projection into [Pk(T)] d , then 

(Q q, Vw) T = (q, Vw) T 

(5.10) = -(V • q,to)T + (q- n,w)ar, 

where we have applied the usual divergence theorem in the second line. Substituting 
(|5.10p into (|5.9|) yields ()5.7|) . This completes the proof of the lemma. □ 

The second equation is concerned with the bilinear form (V<j-v, Lk+iw) for v £ Vh 
and w £ i/ 1 (J7). Using the definition of Lh and the integration by parts, we have for 
any v = {vo, Vf,} £Vh and w £ i? 1 (T) that 

(Vd • v, L h w) T = -(v , V(L h w)) T + (v b • n, L h w) dT 

= (V • v , L h w) T + ((vfc - v ) • n, L h w) 9 T 

= (V • v , w) T + ((v 6 - v ) • n, L h w) aT 

= -(v , Vw) T + (v • n, w) aT + ((vj - v ) • n, L h w) dT 

= -(v , Vw) T + ((v - Vft) • n, io - L h w) dT + (vj • n, w) aT - 

The result can be summarized as follows. 

Lemma 5.3. For any v = {vo, Vf,} £ Vh and w £ iJ x (T) on the element T £ Th, 
one has 

(V d • v, L h w) T = -(v , Vw) T + ((v - Vft) • n, w - L h w) dT 

(5.11) +(v 6 • n, w)st. 

Moreover, if w £ i? 1 (J7), i/ien t/ie sum o/ i5. iij) over all T £Th gives 

b(v, L^w) = -(v 0) V/jw) + ((v - Vfe) • n, w - L h w) dT 

T£T h 

(5.12) +(v b • n, w}dn, 
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where V^iu * s the gradient of w taken element-by- element. 

Let T G 7?,, be any element, and q £ [H 1 (T)] d . Then, we have from the definition 
of Qb that 

||q-n-<9 6 (q-n)||| T = (q • n - Q b (q • n), q • n - Q b (q • n)) dT 
= (q n Q b (q • n), q • n (Q q) • n)) aT 
< ||q - n - Q 6 (q-n)|| aT ||q - n - (Q q) • n|| aT . 

It follows that 

(5.13) q n Q b (q • n)|| ST < ||q • n - (Q Q q) ■ n\\ dT . 

Lemma 5.4. Let q £ [H k+1 ($l)] d and s be any real number such that < s < k. 
Then, we have 

(5-14) M|q-n-Q b (q.n)||! T < h 2 ^ \\q\\ t+1 , 

Ten 

and 

(545) H-QhCli < ph s+1 \\q\\ s+1 . 



Proof. Apply the trace inequality (IA.1I) (see the Appendix) to the righ-hand side 
of (|5.13p to obtain 

(5.16) h T \\q • n - Q b (q • n)||| T < ||q - Q q|lT + h 2 T \\V(q - Oo«l)||r- 

Thus, it follows from Lemma [541 that the estimate (|5 . 14[) holds true. 

To derive (|5.15[) . we have from the definition of the bilinear form s(-, •) in (|3.12l) . 
the definition of Q b , and the inequality (|5.13j) that 

s(q - Qhq,q- Qhq}= p Y MKq- Qoq) • n - (qn- Q 6 (q- n))||| T 

Ten 

<pJ2 ft T (||(q-Qoq) •n||| T + ||q-n-Q 6 (q-n)||| T ) 

Ten 

^pY Ml(q-Qoq) -Hllr- 

TGT h 

Now apply the trace inequality (|A. 1|) to the last term of the above inequality to obtain 

(5.17) s(q - Q h q, q - Q h q) < p Y (||q - Q q||f, + h 2 T \\ V(q - Q Q q)|||) . 

TdT h 

Finally, it follows from the triple-bar norm definition (|4.1j) and the above inequality 
that 

lllq - Q/iq||| 2 = a(q - Qhq, q - Qhq) + s(q - Q h q, q - Qhq) 
= \\®Hq - Qoq)\\ 2 + s(q - Qhq, q - Qhq), 



which, combined with Lemma 15.11 gives rise to the desired estimate (I5.15|) . This 
completes the proof. □ 
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6. Error Analysis. The goal of this section is to derive some optimal order error 
estimates for the WG-MFEM approximation (qh',Uh) obtained from (|3.13[) - (|3.14j) . 
The triple-bar norm as defined in (14.11) will be used to measure the error for the 
vector- valued function q — q^ , while the discrete H 1 norm of (|4.2[) shall be applied to 
u — Uh- Note that the triple-bar norm is essentially a discrete L 2 norm. 

In finite element analysis, it is routine to decompose the error into two components 
in which the first component is the difference between the exact solution and a properly 
defined projection of the exact solution while the second component is the difference 
of the finite element solution with the same projection. In the current application, we 
shall employ the following decomposition: 

q - q h = (q - Q h <i) + (Q h q - q h ), 
u-Uh = (u - L h u) + (L h u - Uh). 

For simplicity, we introduce the following notation: 

eft = {e , e b } := q h - Qhq, eh ■= u h - L h u. 

Our objective is to establish some estimates for e/j in the triple-bar norm and e/, in 
the discrete H 1 norm. For simplicity of notation, we shall denote X^Ter II^HIt by 
\\V h w\\ 2 for any w G Uren H ^ T )- 

6.1. Error equations. Let v = {vo,V(,} G Vh be any finite element function. 
Testing the first equation of (jl.ljl against vo gives 

(aq,v ) + (Vti,v ) = 0. 
Now applying the identity (|5.12j) to the term (Vw,Vo) yields 

a(q, v) - 6(v, L h u) = ^ (( v o ~ v fc) ■ n ; L hU - u) dT + (g, v 6 • n) a n, 

T£T h 

where the Dirichlet boundary condition (jl.2l) was also used. Recall that a s (v) = 
a(-, •) + s(-, •). Adding and subtracting a s (Qhq, v) in the equation above gives 

a s (Qhq, v) - 6(v, L h u) = a(Q h q - q, v) + s(Q h q, v) + (g, v b • n)on 

(6.1) + 7^ ((vq - Vh) • n, L h u - u)ar- 

T£T h 

It follows from the definition of the bilinear form s(-, •) that for any w G [H 1 (il)] d , 
one has 

(6.2) s(w,v)=0, VvgV. 
Thus, we have s(Qhq, v) = s(Qhq — q, v) and 

a(Qhq - q, v) + s(Q h q, v) = a s (Q h q - q, v). 
Substituting the above into (|6.ip yields 

a s {Qh^ v) - b(v,L h u) = a s (Q h q - q, v) + (g,v b ■ n) dn 

(6.3) + ^ ((vp - v&) • n, L h u - u)^. 

Ten 
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Next, we test the second equation of (jl.ip against any w G W/, to obtain 



(V-q,to) = (f,w). 



Substituting (|5.8j) into the above equation yields 



(6.4) 



b{Q h ci, w) = (/, w)- ^ (q • n Q h (q • n), w) dT - 
Ten 



Now, subtracting (|6.3|) from (13.131) gives 



(6.5) a s (e h ,v) - b(v,e h ) = a 5 (q- Q^q, v) - ((v - v b ) • n, L h it - u) aT , 



for all v eVii- Similarly, subtracting (|6.4I) from (13.141) yields 



(6.6) 



6(e h , w)= ^ (q ■ n - Q & (q • n), w) 9T , 
Ten 



for all to€>V/ l . 

The equations (|6.5[) and (|6.6[) constitute governing rules for the error term e/j and 
eh,. These are called error equations. 

6.2. Error estimates. Let us first establish some estimates useful in the forth- 
coming error estimates. 

Lemma 6.1. Let u e H k+2 {Vt) and q e [H k+1 (fl)] d be two smooth functions on 
VI. Then, the following estimates hold true 



(6.7) K( v ° - Vfc ) ' n ' LhU ~ u ^t\ < h s+1 \\u\\ s+2 \\w\l 0<s<k, 



(6.8) ]T |(q-n-Q 6 (q-n),^) aT | < ||q|| s+1 ||V^||, < a < fc, 



/or all v G V and tp 6 IlTgTh ^(T). 

Proof. It follows from the Cauchy-Schwarz inequality, the definition of ||| • |||, and 
the trace inequality (IA.1I) that 



Ten 



Ten 






where we have used the estimate (|5.3I) . This verifies the validity of (|6.7[) . 
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Next, let (p be the average of ip over each element T. It follows from the definition 
of Qb, the estimates (|5.14[) . and (|A.1[) that 

E I fa ' n ~ Qbfa ' n ),<P)dT\ = Y \(<l ,n -Qb(<l-n),<p-0)dT\ 

Ten Ten 

< ( Mq-n-06(q-n)llar) ( £ tf\\v-?\\lr) 
\TeT h / \Ten / 

< ^ +1 ||q|| s+1 ( 2 (^T 2 |l^ - p|£ + IIV^IH) J 

\TGTh / 

<ft s+1 ||q|| s+1 ||V^||, 

which completes the proof of (|6 .8[) . □ 

Lemma 6.2. For any 99 G Wh> there exists at least one v € Vh such that 

(6-9) \b{v,<p)\ = h\\l, h , 

(6-10) |||v||| < \\<p\\ lth . 

Therefore, the following inf-sup condition is satisfied 

(6.11) sup > py\\ uh 

vev fe , llvli^o lll v lll 

/or some positive constant f3 independent of the meshsize h. 

Proof. It follows from the definition of b(-, •) and the discrete divergence that 

& ( v >^) = ( V d ■ v > V) 

= X! (( v 6 • n , t P)aT - (v , V<p)r) 

= ^ (v 6 • n e , [^]} e - Y ( v o> V(p) r , 
eefh Ten 

provided that the jump fipj was taken consistently with the direction set T>h- Now 
by taking v = — Vip on each element T and = h~ 1 Qb(lipJ)n. e on each edge or flat 
face e we arrive at 

b(y,<p) = ]T ( V ^' E ^ 1 <06M,Q6M>e = ||v|lU' 

Ten eee h 

which verifies (|6.9[) . To verify (|6.10[) , we apply the definition of the triple-bar norm 
to the above chosen v to obtain 

|||v||| 2 = (aV h <p, V h <p) + P E h T\\y<f ■ n + h~ l Qh{lv\)n e • n]|| T 

<\\^M\ 2 +Y.K 1 \\Q b M\\l = Ml h , 

ee£h 

which is what (|6.10[) states for. □ 

Lemma 6.3. The bilinear form £>(•,•) is bounded in Vh X Wfe. In other words, 
there exists a constant C such that 

(6.12) \b{v,tp)\ <C|||v||| l^llLh. 
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Proof. We note from the definition of weak divergence that 

= 53 i( V b ■ n,<fi)dT - {v ,V(p) T ) 
Ten 

= 53 ( v b ' n e> M)e - J! ( V 0, V^) T - 

eeffc Ten 
Using Cauchy-Schwarz we obtain 

\b(y,<p)\ < K-ne||e IIMIU+ E INHt H v * 

<( 53 Mn- ^112 ) (E^IIMIIe) +]|v || ||v^|| 

(6-13) < H 5> e K-n e || 2 J) +||v ||j |M|u, 

Let T be an element that takes e as an edge or flat face. Then, using the trace 
inequality (|A.1[) and the inverse inequality (|A.21[) we obtain 

Mv 6 ■ n e ||e < 2/i e ||(v fc - v ) • n e f e + 2h e \\v ■ n e ||^ 

<fte||(V6-V )-n e ||2 + ||vo|||. 

Substituting the above inequality into (|6.13|) yields 

|6(v,v)| < |||v||| \\ip\\ lih , 
which completes the proof of the lemma. □ 

With the help of Lemma lBTTI Lemma l6~2l and Lcmma l6~3l we are now in a position 
to derive some error estimates for the WG-MFEM solution (q_h',Uh) in the norm ||| • ||| 
and || • || respectively. 

Theorem 6.4. Let q;, £ Vh and Uh £ Wh be the solution of the weak Galerkin 
mixed finite element scheme iS. 13\) - (3.14\ l- Assume that the exact solution (q;u) of 
(TOP-fl^P is regular such that u £ H S+2 {VL) and q £ [H s+1 {tt)] d with s £ (0,fc]. 
Then, one has 

(6.14) |||qfc - Qhq||| + \\u h - L h u\\ hh < h s+1 (\\u\\ s+2 + \\q\\ s+ i) . 

Proof. Let 

£(v) := a s (q - Q h q, v) - 53 ((v a - v&) ■ n, L h u - u) dT 

Ten 

be a linear functional on the finite element space Vh- Also, let 

XO) := 53 ( q ■ n ( 5 b ( q ■ n )' W ^ 9T 
Ten 
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be a linear functional on the finite element space W^. The error equations (16.51) and 
(|6.6[) indicate that the pair (e/,; e^) is a solution of the following problem 

(6.15) a s (e h ,v)-b(v,e h )=l(v), Vv 6 V h , 

(6.16) b(e h ,w)=xM, VweW h . 

The bilinear form a s (-, •) is clearly bounded, symmetric and positive definite in Vh 
equipped with the triple-bar norm ||| • ||| . Lemma 16.31 indicates that the bilinear form 
&(•, •) is bounded in Vh x Wh, and Lemma 16.21 proved that the usual inf-sup condition 
of Babuska [3] and Brezzi [5] is satisfied. Thus, we have from the general theory of 
Babuska and Brezzi that 

(6-17) l|efc||| + ||eh||i,h<c(||^|| n + ||xllwi), 

where \\t\\v and ||x||w' are the norm of the corresponding linear functionals. To 
estimate the norms, we use (|5.15j) and (|6.7|) to come up with 

K(v)|<^ +1 (llq||, + i + ll«l| s+ 2)IHI|. 

Thus, we have 

(6-18) \\e\\ Vh <h 8+1 (\\q\\ a+1 +\\u\\ a+2 ), 0<s<fc. 

As to the norm of x, we use the estimate (16. 8[) to obtain 

IxNI </* s+1 ||q|Ui||H|i )/s , 0< S <fc, 

which implies that 

(6-19) \\x\\w' h <h s+1 \\q\\ s+ i, 0<s<k. 

Substituting (16.181) and (|6 . into (|6 . 1 7[) yields the desired error estimate (|6 . 14|) . □ 

7. Error Estimate in L 2 . To obtain an optimal order error estimate for the 
scalar component e% — Uh — LhU in the usual L 2 -norm, we consider a dual problem 
that seeks \1/ and <fi satisfying 

(7.1) a* + V0 = O, inO 

(7.2) V • * = e h , in tt 

(7.3) = 0, on dCl. 

Assume that the usual _ff 2 -regularity is satished for the dual problem; i.e., for any 
£/, G L 2 (Q), there exists a unique solution (^; </>) G [// 1 (il)] d x H 2 {Vt) such that 

(7-4) Ml2 + ||*||i<IM|. 

Lemma 7.1. Let (q;u) be the solution of tl.3\) - [L4\ ). Assume that u e H k+2 (ft), 
q € [H k+1 (fl)] d , and the coefficient junction satisfies a € W 1,ca {T) on each element 
T . Then, one has the following estimates 

(7-5) Mq-Qhq,Q^)|<^' +2 ||q|U+l||*||l, 

(7.6) ]T KQo^-n-QiC^.nJ.LfcU-^arl^/i^Hlfc+all*!!!, 

Ten 

(7.7) Kq-n-Q fc (q-n),L^-0) aT |<^+ 2 ||q|| fc+1 ||0|| 2 . 

TeT h 
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Proof. Let a be the average of a on each element T. It follows from the definition 
of a s (-, •) and the approximation property f|5 . 1 5[) that 

\a s (q-Q h q, Q h ^)\< |(a(q - Q q), Q *)| + |»(q - <2oq, 

= l(q-Qoq, (a-a)Q *)| + |*(q-0oq,0fc*-*)| 

< C^IIVallooHq- Qoq|| ||Q„*|| + |||q - Q h q||| |||* - Q h *||| 

S/^Hqllk+ilMU, 

where ||Va||oo is the L°° norm of Va taken on each element. As to (|7.6[) . we use (|A.1|> 
and (|5. 15|) to obtain 

E |(Q tf •n-Q 6 (*-n) 1 Z h tt-tt} ar | 
Ten 

= \((Qo^-^)-n-(Q b (^-n)-^-n),L h u-u) dT \ 

Ten 

\ 1/2 

- it||| + \\V(L h u - u)\\t) J 



\Ten 



<h k+2 \\u\\ k+2 \\H>\\i- 
Similarly, using (15. 14|) we have 

E Kq-n-Q^q-n): L h ct> - 4>} dT \ 

Ten 

< ^ +1 ||q|U+i ( E (^ 2 H^ - ^IIt + HV(^ - 0)111) ] 
\Ten ) 

<^ +2 ||q|U + i||0|| 2 . 

This completes the proof of the lemma. □ 

We are now in a position to establish an optimal order error estimate for the 
scalar/pressure component eh = it/, — L^u in the usual L 2 -norm. The result is stated 
in the following theorem. 

Theorem 7.2. Let q h e V h and u h e W h be the solution of \3.13)) - ^TT$ . As- 
sume that the exact solution of H1.3\) - [T^ satisfies u G H k+2 (Q) and q G [H k+1 (Q)] d . 
Then, one has the following error estimate 

(7-8) |t M , l -L /lU ||<^+ 2 (|| M || fc+2 + ||q|| fe+1 ). 

Proof. Testing (|7.2I) by e/, = — and then using (|5.8[1 yields 

(e fc , e h ) = (V • *, e h ) 

(7.9) = 6(Q fe *, e h ) + E <* ' n - g b (* • n), e h ) dT . 

Ten 

It follows from (16.51) that 



+ E((<2o*) ' n - Q*>(* ' n), Lftit - u) 9T . 
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Substituting the above equation into (|7.9p yields 

(en, e h ) = a s (e h , Q h ^) - a s (q - Q h q, Q h ^) 

(7.10) + Wo* • n - Q b (* • n), L h u - u) 9T 

Ten 

+ Y (*-n-Q 6 (*-n), e h ) dT . 

Ten 

Using (|7.6[) we obtain 



(7.11) Yl (Qo* ■ n - ■ n )' LhU - u ^t < h k+2 \\u\\ k+2 11*11!. 

Ten 

Using (|6.8[) we arrive at 

(7-12) £ <# ■ n - Q 6 (* • n), e h ) dT < h\\V h e h \\ ||*||i. 

From (17.51) we have 

(7-13) |o,(q-Q h q, Q h *)| < ^ +2 ||q|| fe+ i 

Now substituting (f7TTT]) - (f77T3|) into ([7J0)) we obtain 

(7.14) (e h , e h ) < \a s (e h , Q h *)| + ^'+ 2 (||q|| fe+ i + ||u|| fe+2 ) ||*|| a + /i||V h e h || ||*||i. 
It remains to deal with \a s (eh,Qh^)\ m (|7. 14|) . To this end, we note that 
a s (e h , Q h ^f) = a s (e h , Q h V - *) + a s (e h , 

and 

|a a (e„, (5^* - *)| < |||e h ||| |||Q h ¥ - *||| < ft|||e fc ||| ||*||i. 

Thus, 

(7-15) \a s (e h ,Q h ^)\<\a s (e h ^)\+h\\\e h \\\ 

Since * € [i? 1 (Sl)] d , then we have a s (e^,*) = a(eft,, i.e., the stabilization term 
vanishes when one of the components is sufficiently regular. Now testing (|7.ip against 
e h gives 

a(e fe ,*) = (a*,e ) = -(V0,e o ). 

Furthermore, we have from (|5.12[) (with v = and w — <f>) and the fact that = 
on dil that 

(V<£, e ) = -&(efe, 2^0) + ((e - e&) • n, <f> - L h 4>) dT . 

TeT h 

Using the error equation (I6.6[) to replace the term b(eh,Lh(f>) in above equation we 
obtain 

(V0, e ) = - ^ (q • n Q 6 (q • n), L h cj)) dT + ^ (( e o - e&) • n, - L h cj)) dT 
Ten Ten 

= ( q ' n < 5 fe ( q ' n )' ^ ~ L h4>)dT + ((eo - e 6 ) ■ n, <f> - L h (j>) dT , 
Ten Ten 
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where we have used the fact that (f> G and <f> = on dil in the second line. 

Now using (|7.7p to estimate the first summation in the above equation, and (|6.7[) with 
s = 0, v = e/j, it = to estimate the second summation we obtain 

\(V^e )\<(h k+2 \\ q \\ k+1 +h\\\e h \\\) \\<P\\ 2 . 

Thus, 

|o.(e fcj *)| - |o(e fcl *)| - |(V0,eo)| < (/i fc+2 ||q|| fe+ i + ft|||e h |||) ||0|| 2 . 
Substituting the above into (17. 15)) yields 

(7.16) |o a (e h> Q J ,*)|<(ft fc + 2 ||q||j k +i + /i|||e h |||) \\cb\\ 2 + h\\\e h \\\ \\*\\i- 
Now substituting f|7. 1 6[) into f|7. 14[) gives 

(7.17) (e h , e h ) <(^' +2 ||q|U+i + h\\\e h \\\) ||0[[ a + h\\\e h \\\ 

+/ i fc+2 (||q|U + i + ||u|| fe+2 )||*|| 1 + / l ||V fe e h || II^Hx. 

Observe that the iJ 2 -regularity of the dual problem implies the following estimate 

ll* + ll*lli<IM|. 

Substituting the above into (|7.17|) and then dividing both sides by \\eh\\ we obtain 

INI < ft(||e fc || + ||V h e h ||) + h k+2 (||q|U +1 + HuHfc+a) , 

which, together with (|6.14j) . implies the desired L 2 error estimate (|7.8|l . This com- 
pletes the proof. □ 

Appendix A. The goal of this Appendix is to establish some fundamental esti- 
mates useful in the error estimate for WG-MFEM. First, we derive a trace inequality 
for functions defined on the finite element partition Th with properties specified in 
Section [H 

Lemma A.l (Trace Inequality). Let Th be a partition of the domain VL into poly- 
gons in 2D or polyhedra in 3D. Assume that the partition Th satisfies the assumptions 
(Al), (A2), and (A3) as specified in Section^ Then, there exists a constant C such 
that for any T £ Th and edge/face e £ dT , we have 

(A.i) \\9\\ 2 e <c{h^\\e\\ 2 T + h T \\ve\\ 2 T ), 

where 9 6 H 1 (T) is any function. 

Proof. We shall provide a proof for the case of e being a flat face and 9 € C 1 (T). 
To this end, let the flat face e be represented by the following parametric equation: 

(A.2) x e = <p(Z,n) :={Mt,v),Mt>ri),Mt>ri)) 

for (£, rf) S D C K 2 . By Assumption A3, there exists a pyramid P(e, T, A e ) with apex 
\,X2,x^) contained in the element T. This pyramid has the following 
parametric representation: 

(A.3) x(i,£, 77) = (1 77) +tx* 
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for (t,^,rj) G [0, 1] x D. For any given x e G e, the line segment joining x e and the 
apex x + can be represented by 

x(t) = x e + £(x* - x e ). 
From the fundamental theorem of Calculus, we have 

6> 2 (x e ) - 9 2 {x(t)) = - ( d T e 2 {x e +TUj)dT, w = x*-x e . 
Jo 

The above can be further rewritten as 

2 (x e ) -9 2 (x(t)) = -2 f 9{V6-uj)dT. 

Jo 

It follows from the Cauchy-Schwarz inequality that for any t £ [0, ^] we have 

i i 

9 2 {x e ) < 6 2 (x(t)) + [ 2 6 2 \uj\dT + hr I 2 |V6»| 2 |w|dr. 

Jo Jo 

Now we integrate the above inequality over the flat face e by using the parametric 
equation (|A.2j) . yielding 



|2 (x e )|0 c x <£„|d£d77 < / 9 2 {x{t))\^ x <t> n \d£d-n 

D J D 

(A.4) + h T 1 I I 2 \uj\\^ x cf> v \d£,dT]dT 

Jo Jd 

1 

(A.5) +h T \ [ |V(9| 2 M|0 C x tnldtdrjdT. 

Jo Jd 

Observe that the integral of a function over the following prismatoid 
Pi :={x{t,£,Ti) : (*, ^, 77) G [0,1/2] x D} 

is given by 



/(x)dx = / / /(x(r^,7 ? ))J(T,C^)^^ ) 

Pi JO Jl5 

2 

where J(t, f, 77) = (1 — r) 2 |(^ x • w| is the Jacobian from the coordinate change. 
The vector cj>^ x (f> n is normal to the face e, and uj = x* — x e is a vector from the base 
point x e to the apex x*. The angle assumption (see Assumption A3 of Section [3]) for 
the prism P(e,T, A e ) indicates that the Jacobian satisfies the following relation 

(A.6) J(r,f,i?)>j^x^| H, re [0,1/2] 

for some fixed no G (0,1). Observe that the left-hand side of (|A.4j) is the surface 
integral of 9 2 over the face e. Thus, substituting (|A.6j) into (|A.4|) and (|A.5|) yields 

? 2 de< / 6 2 (x{t))\<f> ( x (jy^dtdr) 
Jd 

+4fi 1 h T 1 / 9 2 dx + A^ 1 h T I \V9\ 2 dx. 
JPI Jpl 
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Now we integrate the above with respect to t in the interval [0, |] to obtain 
(A.7) i [e 2 de< [ 2 f 6> 2 (x(t))|<^ x ^dtdrjdt 

2 J e Jo J.D 

+2^ 1 h T 1 I 8 2 dx + 2 f iQ 1 h T [ |V6»| 2 dx. 
Jpi J Pi 

2" 2 

Again, by substituting (|A.6[) into the right-hand side of (|A.7|) we arrive at 

(A.8) i / 6 2 de < 4/ioVr 1 / 2 / 9 2 (x(t))J(t,(,r 1 )d^dr]dt 

+2^ 1 /i T 1 / 6 2 dx + 2^ 1 h T I |V6>| 2 dx. 

2 2 

The integral on the right-hand side of (IA.8I) is the integral of 8 2 on the prismatoid 
Pi . It can be seen from the Assumption A3 that 

< u^hTf} 

for some positive constant a*. Therefore, we have from the above estimate that 
/ Q 2 de<CK T 1 [ 6 2 dx + Ch T [ |V(9| 2 dx, 

Je J Pi JPi 

which completes the proof of the Lemma. □ 

Next, we would like to establish an estimate for the L 2 norm of polynomial func- 
tions by their L 2 norm on a subdomain. To this end, we first derive a result of similar 
nature for general functions in H . 

Lemma A. 2. Let K C M. d be convex and v G H l (K). Then, 

(A.9) \H 2 K < j^-hf s + -^-\\^f K , 

where 5 is the diameter of K and S is any measurable subset of K. 

Proof. Since C (K) is dense in H 1 ^), it is sufficient to establish (IA.9I) for 
v G C X {K). For any x, y £ K, we have 

/■l x -y| v-x 
v 2 (x) = v 2 (y) - / d r v 2 (x + ruj)dr, uj = f r . 

Jo |y - x l 

From the usual chain rule and the Cauchy-Schwarz inequality we obtain 

/■l x -y| 

v 2 (x) = v 2 (y) — 2 / vd r v(x + rui)dr 
Jo 

(A.10) <v 2 (y)+e5- 1 [ l v 2 dr + e^S /' \d r v\ 2 dr, 



where e > is any constant. Let 



w 2 (x), x G K 

0, x 4 K 
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\d r v\ 2 (x), x e K 

0, xiK. 



and 

W(x) -- 

Then, the inequality (|A.10[) can be rewritten as 

pOO />OC 

v 2 (x) < v 2 (y) + eS^ 1 V {x + ruj)dr + (7 1 5 j W(x + ru)dr. 
Jo Jo 

Integrating the above inequality with respect to y in S yields 

(A.ll) \S\v 2 {x) < f v 2 dS 
Js 

+ ( eS l V{x + ruj)dr + e~ 1 6 W{x + ru)dr J dy. 

J\x-y\<8 \ Jo Jo J 

It is not hard to see that 

f5 



p p OO p OO p po 

/ / V(x + ruj)drdy = / / / 

J |x— y|<5 Jo Jo J \u)\ — l J 



V{x + ruj)p d 1 dp dco dr 



S d f f 

— / V(x + rw)doj dr 
d J J\ u \=i 



(A.12) | x _ y |i-V(y)rfy. 
Analogously, we have 

(A.13) f H W{x + ru)drdy= 5 - [ \x - y^dMy^dy. 

J\x-y\<6J0 JK 

Substituting ([A~T2|) and (|A~T3)) into (fA~TTj) yields 

\S\v 2 (x)< f v 2 dS+^ f | x - y | 1 -V(y)dy + ^ f |x - yM V«(y)| 2 dy. 

Js « JK ea Jk 

Now integrating both sides with respect to x in K gives 

v 2 dK<\K\ [ v 2 dS+— [ v 2 dK+^- [ \\7v\ 2 dK, 
k Js d Jk Jk 

which yields the desired estimate (|A.9[) by setting e = . □ 



Consider a case of Lemma IA.2I in which the convex domain if is a shape regular 
d-simplex. Denote by hx the diameter of K . The shape regularity implies that 

1. the measure of K is proportional to , 

2. there exists an inscribed ball Bk C K with diameter proportional to \ik- 
Now let 5* be a ball inside of K with radius r$ > s*hx- Then, there exists a fixed 
constant such that 

(A.14) \K\ < k*\S\. 
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Apply (|A.14|) in (|A.9[) and notice that r»|5| > h d K and (5 = Thus, 

(A.15) | |t ,||^<2«.|| w ||| + ^^:|[Vt,||^ 

For simplicity of notation, we shall rewrite (|A.15[) in the following form 

(A.16) IMI^<aolMls + ai/4l|V*. 

If v is infinitely smooth, then a recursive use of the estimate (|A.16I) yields the following 
result 

n 

(A.17) || V ||^<^ a3 ^||V^|| 2 s + a„ +1 ^+ 2 ||V"+ 1 H|^. 

3=0 

In particular, if v is a polynomial of degree n, then 

n 

(A.18) [[*<E^I|VMI- 

3=0 

The standard inverse inequality implies that 

l|V^|| s < h-J\\v\\ s . 
Substituting the above into ([A. 181) gives 

(A.19) IMIi<IMII- 

The result is summarized as follows. 

Lemma A. 3 (Domain Inverse Inequality). Let K C IR d be a d-simplex which has 
diameter hx and is shape regular. Assume that S is a ball in K with diameter r$ 
proportional to hx; i.e., r$ > ^*hx with a fixed <r* > 0. Then, there exists a constant 
C = C(<j*,n) such that 

(A.20) \M 2 K<C(<;*,n)\\vf s 

for any polynomial v of degree no more than n. 

The usual inverse inequality in finite clement analysis also holds true for piecewise 
polynomials defined on the finite element partition Tk provided that it satisfies the 
assumptions A1-A4 of Section [31 

Lemma A. 4 (Inverse Inequality). LetTh be a finite element partition of Q con- 
sisting of polygons or polyhedra. Assume that 7~h satisfies all the assumptions A1-A4 
as specified in Section^ Then, there exists a constant C — C(n) such that 

(A.21) \\V<p\\ T < C(n)/^ 1 IMh VTeT h 

for any piecewise polynomial ip of degree n on 7h- 

Proof. The proof is merely a combination of Lemma lA.31 and the standard inverse 
inequality on d-simplices. To this end, for any T £ Th, let S(T) be the circumscribed 
simplex that is shape regular. It follows from the standard inverse inequality that 

||V^||r < ||V^|| s( t) < Ch^Msix). 
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Then we use the estimate (|A.20|) . with K = S(T), to obtain 

\\V<p\\T<Chj. 1 \\(p\\s<Chj. l \\<p\\ T , 

where S is a ball inside of T with a diameter proportional to hx- This completes the 
proof of the lemma. □ 
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